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As a minimal fermionic model with kinetic frustration, we study a system of spinless fermions in 
the lowest band of a triangular lattice with long-range repulsion. We find that the combination of 
interactions and kinetic frustration leads to spontaneous symmetry breaking in various ways. Time- 
reversal symmetry can be broken by two types of loop current patterns, a chiral one and one that 
breaks the translational lattice symmetry. Moreover, the translational symmetry can also be broken 
by a density wave forming a kagome pattern or by a Peierls-type trimerization characterized by 
enhanced correlations among the sites of certain triangular plaquettes (giving a plaquette-centered 
density wave). We map out the phase diagram as it results from leading order Ginzburg-Landau 
mean-field theory. Several experimental realizations of the type of system under study are possible 
with ultracold atoms in optical lattices. 



Geometric frustration in classical and quantum many- 
body systems is a source of intriguing phenomena like 
extensive ground-state entropies, topological order, and 
exotic emergent low-energy physics [T]-[5]. It is naturally 
encountered in systems of antiferromagnetically coupled 
localized magnetic moments, arranged in a non-bipartite 
(e.g. triangular) lattice geometry that prohibits the fa- 
vored antiparallel orientation between all pairs of neigh- 
boring moments. Recently, geometric frustration has also 
been induced to the kinetics of a system of ultracold 
bosonic atoms on a triangular lattice by dynamically in- 
verting the sign of the tunneling matrix elements [6l [7] . 
For weak interaction the system shows spontaneous time- 
reversal (TR) symmetry breaking, whereas strong inter- 
actions are conjectured to lead to spin-liquid-like quan- 
tum disordered behaviour. 

Here we investigate a minimal fermionic model with 
kinetic frustration: spinless fermions on a triangular lat- 
tice with off-site repulsion. It can be realized, e.g., with 
ultracold dipolar atoms or molecules in an optical lat- 
tice [SHTO]. In this system kinetic frustration appears 
naturally as a consequence of fermi statistics. Namely, 
for filling well above one half the system is governed by 
the low-energy states of the fermionic holes whose kinet- 
ics is determined by the sign-inverted tunneling matrix 
elements. We focus on the weakly interacting regime 
and find that, depending on the filling, two types of loop 
currents can spontaneously emerge at experimentally ac- 
cessible temperatures and break TR symmetry. One of 
them breaks the inversion symmetry of the lattice and is 
chiral, while the other one breaks the translational lattice 
symmetry. TR symmetry is broken neither as a conse- 
quence of coupling the kinetics to further degrees of free- 
dom (such as spin or sublattice orbital freedom in the unit 
cell) nor in a process involving (quasi) long-range order in 
a continuous degree of freedom (i.e. bose condensation). 




FIG. 1. (color online) a) and b) Phase diagram resulting 
from low-order Landau expansion of the mean-field free en- 
ergy, with respect to interaction strength [//t, temperature 
T/t, and filling v. Solid (dashed) lines indicate second (first) 
order transitions; T/t = in a) and [//t = 2 in b). c) Frus- 
trated triangular lattice, with solid (dashed) lines represent- 
ing negative (positive) hopping matrix elements, d) Pattern 
of staggered chrial plaquette currents (SC). e) Translational 
symmetry breaking by modulated currents (MC, arrows), a 
density-wave (DW) defining a kagome lattice (larger dot rep- 
resents increased density), and Peierls-type trimerization [PT, 
triangular plaquettes with increased (lowered) correlations are 
indicated by dark (light) green shading.]. 



Typically at least one of these two themes in encountered 
when TR symmtry is broken in fermionic systems; exam- 
ples range from ferromagnetic metals to Mott-insulators 
with spin or orbital order, but comprise also exotic states 
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like chiral spin liquids (not featuring bose condensation) 
pT] OT p ^ ip superconductors (not requiring spin) [12] - 
fT5] . We also observe that the approximate nesting prop- 
erty of the fermi surface near hole filling of 1/4 gives rise 
to a rich spectrum of three different types of instabili- 
ties that break the translational symmetry: One of them 
leads to spatially modulated currents (MC) and has been 
mentioned already. The other two give rise to a density 
wave (DW) and to a Peierls-type trimerization (PT, a 
plaquette-centered density wave). 

Model Our system is governed by the Hamiltonian 
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'=ribe, 
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Here denotes the annihilation operator for a fermion at 
the lattice site located at r, and the vectors connect 
neighboring sites (see Fig. [l]^). Repulsive interactions 
V{t) between distant sites can be achieved in various 
ways, e.g. by using dipolar atoms or molecules that are 
polarized perpendicular to the lattice plane [H fT6H2Q] or 
as a superexchange process between neighboring sites in a 
mixed Mott insulator of fermions and bosons [2TH24] . For 
simplicity, in the following we assume isotropic nearest- 
neighbor interactions of strength U = y(±ei) > 0. The 
tunneling matrix elements along the three lattice direc- 
tions ei shall be given by t\ = —t and t2 = ts = t with 
t > (see Fig. [l]^), giving a 7r-fiux through each plaque- 
tte. This sign configuration corresponds simply to study- 
ing holes instead of particles (which makes a difference 
in the non-bipartite triangular lattice) in a momentum 
shifted reference frame [25]. Another option would be 
to control the tunneling parameters via lattice shaking 
[6l [7] . The system of N particles on M sites (with filling 
ly = N/M) is characterized, moreover, by the tempera- 
ture T. 

The tunneling coefficients are such that the kinetic en- 
ergy cannot be minimized between all neighboring sites 
at the same time, since negative (positive) ti favor a 
single particle state with a relative phase of tt (0) be- 
tween two sites; the system is kinetically frustrated [6]. 
As a consequence, the single-particle dispersion rela- 
tion £{k) = cos(ei • k) possesses two minima (see 
Fig. [2|i). It is possible, however, to maximize the kinetic 
energy on every pair of neighboring sites, so the system 
is sensitive to the frustration for low filling only. 

The double-well structure of the dispersion relation 
supports various interaction-driven instabilities. One 
type of instability results from interaction favoring par- 
ticles to be close by in (quasi) momentum: The potential 
energy between two fermions with sharp momenta p and 
p+gr reads V{0) — V{q) = V{q)^ where the Fourier trans- 
form of the interaction V{q) = U cos(ei • q) de- 
cays with \q\ (see Fig.|2]3), and the relevant momentum- 
dependent second term (the exchange term) obtains a 
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FIG. 2. (color online) a) Single particle spectrum with two 
minima as a consequence of frustration, b) Momentum repre- 
sentation of the nearest neighbor interaction of strength U . c) 
Momentum space density distribution for the staggered cur- 
rent phase (see Fig.[l]l for the real-space signature), d) First 
Brillouin zone with nested Fermi surface (dashed 'bowtie'- 
shaped line) and nesting vectors (white arrows); black line 
indicates reduced Brillouin zone. 



minus sign from fermi statistics. For low filling this can, 
despite the kinetic energy cost, favor an imbalanced oc- 
cupation of the two wells (see Fig. ^) corresponding to 
staggered loop currents (SC) (see Fig. [l]i). 

Another source of instability appears near quarter fill- 
ing V = 1/4: where both minima are filled up completely, 
and the fermi surface is shaped like a bow-tie (Fig. [2]i, 
dashed line) and becomes approximately nested. Namely, 

— ^ sf — £k-\-Q with fermi energy ef and the same 
nesting vector Q for a significant fraction of momenta k 
near the fermi surface, making the kinetic energy cost for 
spatial modulation described by Q low [26^ . The nesting 
occurs for the vectors Qj (Fig. [2]i, white arrows) labeled 
j = 1, 2, 3 such that Ci ■ Qj = (1 — Sij)7T (since 2Qj is a 
reciprocal lattice vector Qj and —Qj are equivalent). We 
find three different instabilities that break translational 
symmetry (see Fig. [l^) to be described below. 

Method We use mean-field theory. In momentum rep- 
resentation, dk = M~^/^ e~*'^ '^c-p, the interaction 
consists of terms V {q)d]^_^^dp_qdpdk that we approx- 
imate like d'ldldpdq ^ (,kqo\dp + iipd]^dq - ^kq^ip - 
ikpo\dq - ^iqd]^dq + ^kp^iq whllc asklug for self consis- 
tency ^kp = (^fc^p)- Allowing for spatial modulations 
described by the nesting vectors, we retain averages of 
the form {d]^_Q,dk) with j = 0,1,2,3 where Qo = 0- 
Accordingly, we find single-particle correlations 
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FIG. 3. (color online) Phase boundary of parameter region 
where staggered currents are expected (from low-order Lan- 
dau expansion of the mean-field free energy) Solid (dashed) 
lines indicate second (first) order transitions, a) U/t = 
1.5, 1.6, 1.8, 2, from small to large regions, b) T/t = (blue, 
crossing ^/-axis) and T/t = 0.2 (red, not crossing ^-axis). 



where Cj 



i{p-Qj)-ei 



i^i-Qj^p) (summing 



over one Brillouin zone). Defining eo = 0, the lo- 
cal density rir = (c|.c-r) is related to the four aver- 
ages coj = pj. The current from r to r + e^, namely 
Ji{r) = 2ti Im {cl_^^.Cr) /h, is described by ImQ^oj = lij 
(using the equivlence of Qj and —Qj). We also define 
Re Ci^oj = Rij for the real (non-current-generating) part 
of the correlations (cj.^g .c-r). The terms po = N/M and 
RiO are non-zero already without symmetry breaking. 

Staggered currents Well below quarter filling, nesting 
is irrelevant and only terms with j = are important. 
The free energy reads 
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with Uk = {o]^^k) and dispersion relation uj^ = 
2 Ei=i [iti + URio) cos{k ■ a) + Ulio sm{k ■ e,)] • If 
F is minimal for finite momentum imbalances 7^0 = 
Ep sin(p - ei)np^ loop currents will emerge. 
In order to learn when this happens, we apply pertur- 
bation theory within the imaginary time formalism and 
Landau expand the free energy in powers of the /^o- 
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where = (^105^205^30) (with power n and transpo- 
sition T) and where M and K are symmetric matrices. 
A second-order phase transition into a phase with finite 
averages Ii = (—1,1,1)/ occurs when m2, the lowest 
eigenvalue of M, becomes negative, while higher-order 
contributions to F remain positive. The resulting phase 
boundaries are plotted as solid lines in Fig. |3] If the 
quart ic term, with Kij = Sij\r]\-\- {1 — Sij)(^ becomes neg- 
ative for 2( < —\r]\ on the line defined by m2 = 0, the 
phase transition becomes first-order (indicated by dashed 
lines in Fig. |3|. Minimising F exactly at a few points in 
parameter space shows that the shift in the first-order 
phase boundary away from 1712 = 0, which enlarges the 
symmetry broken region, is very small (it is not taken 
into account in Fig. [3| . 



The current pattern that we find in the symmtry- 
broken phase at low densities is defined by Ji{r) = 
—2tiIio = J = ±|J|, giving a staggered pattern of cir- 
cular plaquette currents as shown in Fig. [T]i. The state 
is two-fold degenerate and chiral in the sense that both 
TR and lattice inversion transform one solution into the 
other [11]. Experimentally this state can be identified 
by the momentum imbalance visible in time-of-flight ab- 
sorption images. 

The shape of the phase boundary (Fig. [3| can be un- 
derstood as follows. For lower densities, the slope of the 
single-particle spectrum at the Fermi surface is flatter 
(see Fig. |2^), which explains why the critical interaction 
strength can be reduced when the filling is lowered. At 
high densities, moreover, the double-well structure of the 
dispersion relation (i.e. kinetic frustration) is not rele- 
vant anymore. If the filling is too low, however, a finite 
temperature T may overcome the energy scale Uv of the 
interaction and destroy the currents altogether. Gen- 
erally, increasing the temperature counteracts the clus- 
tering effect from momentum-space attraction, making 
a momentum space imbalance less favorable. We have 
performed exact diagonalizations for a small system (3 
particles on 45 sites) at T = 0, and find an increased ad- 
mixture of kinetic energy eigenstates with large \Iio\ for 
/7/t > 1, in support of the mean- field findings. 

Spatial modulation Near quarter filling, where nest- 
ing becomes relevant, we have to take into account the 
additional possibility of spatial modulation, leading to 
a four-site unit cell. The mean-field Hamiltonian pos- 
sesses four bands in a reduced Brillouin zone that has 
the same shape but only half the extent of the original 
one (Figj2]i). Broken translational symmetry is indicated 
by finite averages pj, R^j, or lij (with j ^ 0) that open 
a gap between the lowest bands [27 . Expanding the free 
energy also with respect to those averages, in the leading 
(quadratic) order we find 
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with vectors Dj = (pj^Rjj)^ Rj = {Rkj^Rij) and Jj = 
{Ik,j^ Iij) such that j ^ k ^ I under cyclic permutation, 
symmetric 2x2 matrices A, B and C, and coefficient 
A > 0. The latter implies that Ijj = is favored. The 
remaining averages are coupled pairwise, with different 
nesting directions j uncoupled in quadratic order. 

A DW instability is indicated when the lowest eigen- 
value of A becomes negative. This defines the solid 
boundaries plotted in Fig. [4j In third order we find 
a term Tpip2p2> that couples the density modulations 
along the three directions j. It favors the four- fold de- 
generate kagome pattern shown in Fig. [l^ with reduced 
correlations among the low-density sites due to the cou- 
pling to Rjj. The pattern with high and low occupa- 
tions interchanged is disfavored (an obvious consequence 
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FIG. 4. (color online) Different types of instabilities that 
break the translational lattice symmetry are found inside 
the plotted boundaries. Solid, dotted, and dashed lines 
stand for instabilities leading to density- wave order, Peierls- 
type trimerization (i.e. plaquette-centered density-wave or- 
der), and a pattern of modulated currents, respectively. The 
transition to the density- wave order has to be first order, a) 
U/t = 0.5 (inner curve), 1,1.5,2 (outer curves), b) T/t = 
(lower curves) and 0.1 (upper curves). 

of nearest-neighbor repulsion on the triangular lattice). 
The third-order term also renders the transition first or- 
der and shifts the phase boundary slightly (not taken into 
account in Fig. |4|. Exact diagonalization studies with 4 
particles on 16 sites at T = reveal kagome DW patterns 
with a slight anisotropy (probably a finite size effect). 

If the lowest eigenvalue of B turns negative (Fig. |4j 
dotted lines), an instability towards a Peierls-type mod- 
ulation of nearest-neighbor correlations is found: Choos- 
ing spontaneously an embedded kagome pattern, the cor- 
relations {cl,Cr) on the A (v) triangular plaquettes of 
that pattern are increased (lowered); or vice versa (see 
Fig. [l^). This eight-fold degenerate Peierls trimerization 
breaks the inversion symmetry of the lattice and corre- 
sponds to a plaquette-centered DW. 

Finally, if C acquires a negative eigenvalue (Fig. |4j 
dashed lines) an instability towards the formation of MC 
is indicated. One finds the current carrying bonds de- 
fine a kagome lattice with circular currents of equal ori- 
entation around both types of triangular plaquettes (A 
and v) of opposite orientation around the hexag- 
onal plaquettes (such that the finite plaquette fluxes of 
the mean-field Hamiltonian vanish when averaged over 
the unit cell). The eight-fold degeneracy of this pattern 
results from spontaneous breaking of both translational 
and TR symmetry. 

The structure of the stability boundaries shown in 
Fig. [4] explains itself in an intuitive way. According to 
the nesting condition, interaction-induced spontaneous 
translational symmetry breaking occurs for filling near 
1/4, in a finite interval of filling factors that grows with 
increasing interaction. An interesting effect is that for a 
given interaction strength U/t this interval has its max- 



imum extent at a non-zero temperature. The growth 
of the symmetry broken parameter region when a small 
temperature is switched on, might be explained by the 
fact that smoothening the fermi edge increases the num- 
ber of participating momentum pairs for which the nest- 
ing condition is approximately fulfilled. However, for 
larger temperatures the order is eventually destroyed. 

As visible in Figure [4] there is a large region in pa- 
rameter space, where according to second Landau order 
all three nesting-driven types of symmetry breaking are 
favorable. As long as the spontaneous kagome patterns 
for the three match, they can in principle co-exist; their 
structures do not mutually exclude each other. This sug- 
gests that (at least as long as the order parameters are 
small such that low-order Landau terms dominate) the 
system might take advantage of realizing two or even all 
three of them at the same time, leading to a degeneracy 
of up to sixteen. However, generally it will depend on 
the energetics described by higher-order Landau terms 
that contain the coupling between the different orders, 
whether mean-field theory predicts such co-existences. 
Shifting and/or merging of the phase boundaries due to 
higher-order terms can be expected. Beyond mean-field 
theory this question can be answered by an experiment 
with ultracold atoms. 

All three translational symmetry breaking orders (DW, 
MC, and PT) open up a gap between the lowest two 
bands of the mean-field dispersion relation in the reduced 
Brillouin zone. However, the momentum-space structure 
of the coupling that induces the gap depends on the type 
of order. Thus the three orders can be distinguished by 
the momentum distribution as it can be measured by 
time-of-flight absorption imaging. 

We summarize our findings in Fig. [l] with subfigure a) 
and b) showing the combined phase diagram. For suffi- 
ciently large interactions there is a small parameter re- 
gion where both the SC and the MC overlap. We expect 
both types of order to repel each other and it appears 
likely that higher order terms will remove this overlap re- 
gion and predict a direct transition between both phases. 

We conclude that kinetic frustration can give rise to 
rich physics in fermionic lattice models already for weak 
interactions. As an example we have considered the min- 
imal model of spinless fermions on a triangular lattice for 
which we predict the emergence of chiral or modulated 
loop currents as well as Peierls-type trimerisation and 
density- wave instabilities. Our findings can be probed 
experimentally with ultracold fermionic atoms at acces- 
sible temperatures. Such an experiment would comple- 
ment the recent observation of time-reversal symmetry 
breaking as a consequence of kinetic frustration in a sys- 
tem of bosonic atoms [7^ . 
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